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Abstract. - The energy of an elastic manifold in a random landscape at T = is shown 
numerically to obey a probability distribution that depends on size of the box it is put into. If 
the extent of the spatial fluctuations of the manifold is much less than that of the system, a cross- 
over takes place to the Gumbel-distribution of extreme statistics. If they are comparable, the 
distributions have non-Gaussian, stretched exponential tails. The low-energy and high-energy 
stretching exponents are roughly independent of the internal dimension and the fluctuation 
degrees of freedom. 



The statistical mechanics of elastic objects or manifolds changes in the presence of disorder, 
since temperature becomes irrelevant as a scaling variable. The statistical properties are 
determined by the competition of elasticity and randomness. Examples of systems where this 
happens are domain walls (DW) in random magnets and flux lines in superconductors ^ 
121121 • At the zero temperature "fixed point" in the renormalization group language, the 
geometry of the object becomes critical. The fluctuations are self-affine below the upper 
critical dimension of the Hamiltonian. This upper critical dimension may or may not have 
a finite value depending on the dimensionality of the possible fluctuations. Examples are 
provided by the one-dimensional directed polymer (DP) in n > 1 fluctuation dimensions, 
which has a low-temperature phase for which n c = oo 2 . Another case is given by D 
internal-dimensional random manifolds in d = (D + n) dimensions, where n — 1, which are 
known to have the upper critical dimension D c = 4 [5]. 

The criticality of elastic manifolds is manifested by the distribution of energies, P{E), 
which is not a Gaussian. This was shown by extensive simulations of directed polymers in 
d = (D + n) = (1 + 1) and (1+2) dimensions (d is the total, embedding space dimension of a 
system) 0] . In these simulations a DP was let to minimize its energy by keeping one end fixed, 
and letting the other one wander freely. The outcome is that for energies smaller than the 
average, EL <C (E), the distribution P{E) ~ exp[— i.e., stretched exponential. Note 
that the distribution is normalized in such a way that the average of E is zero and its variance 
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or standard deviation equals unity. Similarly for E + 3> {E) one obtains a stretching exponent 
?7_l_. These stretching exponents are below the Gaussian 7/ = 2 value for small energies and 
above it for large ones, and their numerical values are rather insensitive to whether n = 1 or 
2. 

What has not been addressed so far is how exactly the scaling function, P(E), depends on 
the set of boundary conditions imposed on the manifold and the dimensions n and D. One 
reason is that computing P{E) is naturally a rather formidable task. One rare analogy is 
provided by the (l+l)-dimcnsional Kardar-Parisi-Zhang-equation (KPZ) [Jj], in the form of 
the asymmetric exclusion process (ASEP). This system is governed by the "strong coupling 
fixed point" of the KPZ interface growth. Approaches based on a generalized free energy by 
Derrida and Lebowitz [Sj and on mapping the problem to random matrix theory (by Prahofer 
and Spohn [7]) have succeeded in deriving analytically the non-Gaussian scaling functions. 
This is true for the velocity fluctuations or height fluctuations in the two respective cases but 
only in the steady-state of the growth, that is when the correlation length exactly equals the 
system size. In the KPZ case in (l+l)-dimensions one can also resort to scaling arguments 
for the extremal properties of the interface statistics, predicting rf- =3/2 [Hj. 

In this Letter we show that the problem of the shape of the probability distribution P{E) 
is actually continuously dependent on a natural boundary condition, namely the transverse 
height of the system. It is also, on the other hand, rather independent of the dimensionality 
The continuum of exponents 77+, ?7_ for P(E) results from the cross-over from a Gaussian 
distribution (when the height of the system is limited to below the natural roughness or 
geometrical fluctuations and hence reflects Poissonian statistics) to the Gumbel function of 
extreme statistics [51ll(Jj when the transverse height of the system is let to increase to infinity. 
It turns out to be so that the variation of the rj-, 77+ is related to so-called penultimate 
extremal distributions, the scaling of an extremal quantity from samples with a finite number 
of independent variables. 

There remain several open questions that we return to in the summary. Recall that 
the energy of a DP relates to the first arrival time of a KPZ interface to a fixed height (a 
related problem but not quite the same is the maximal height of an interface with a given 
average height Hence varying the transverse height thus equals moving between the 

stationary state and an initial transient. The former can be accessed via the large deviation 
and random matrix formalisms [H|[7|. In our numerical simulations we have varied both the 
internal dimension of the system D, which gets values from one to three, and the fluctuation 
degree of freedom n, which gets the same values. We are for algorithmic reasons restricted to 
cases, where whenever D > 1 n has to have a value of unity, and vice versa. 

The geometric fluctuations of clastic manifolds can be measured by the two-point cor- 
relation function or by mean-square fluctuations. These define the roughness exponent £, 
e.g. as w 2 = ([z(x) — z(x)] 2 ) ~ L 2 ^, where L is the linear size of the system, z(x) is the 
manifold's height and x is the D dimensional internal coordinate of the manifold. For the 
ground state or free energy fluctuations around the disorder-averaged mean one has that 
AE = ((E - {E)) 2 ) 1/2 ~ L e , where 9 = 2( + D - 2 [T2]. The ground state energy is a 
quenched random variable and it behaves as E ~ AL D + BL . Thus normal equilibrium 
behavior is modified, and the (free) energy develops non-analytic corrections. In a renor- 
malization picture the energy fluctuation exponent 9 arises naturally. We concentrate on 
random-bond disorder, which means that the random-potential in the energy Hamiltonian is 
delta-point correlated, although the results should be applicable to other types of short-range 
correlated randomness, too. At low temperatures in (1 + 1) dimensions, due to the equivalence 
of the directed polymer in random media |5I2| to the KPZ equation, the exact roughness expo- 
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Fig. 1 - A sketch of (D + n) = (2 + 1) dimensional manifold in a system, whose lateral size, perpen- 
dicular to the average normal of the surface, is L x = L y = L and transverse height, parallel to the 
average surface normal, is L z . The width, root-mean-square fluctuation, of the surface is w. When 
the manifold is pinned in the center of the system, the "cone" geometry is used. In that case the area 
where the surface is not allowed to pass through is depicted with gray. 

nent is ( = 2/3, with RB disorder |5I2I13| . In higher D dimensions functional renormalization 
group calculations give £ ~ 0.208(4 — D) The (1 + n) dimensional directed polymers with 
n > 1 have roughness values ( — 0.62 and £ = 0.59, when n = 2 and n = 3, respectively (see 
e.g. PI). 

For numerical calculations, when D > 1, the optimization of the Hamiltonian of a man- 
ifold in random bond environment (which is equivalent with domain wall in random-bond 
Ising magnet) maps into the minimum cut-maximum flow problem of combinatorial optimiza- 
tion |14U15| . The simulations of D = 1, n > 1 directed-polymers are done using the standard 
transfer-matrix (TM) method ^B] . Here we impose two different type of boundary conditions 
in the transverse height direction of the system, i.e., in the direction of the average normal of 
the D dimensional manifold. The manifold is fixed in a single valley by pinning it at certain 
coordinate, or let to float freely, so that the mean height can vary in a box defined by the 
system height L z . For DP's the pinning follows in the TM simulations from having one end 
fixed and the other one free. For D > 1 manifolds we construct a new way to create a "one 
valley" case since the same trick as for DP's to fix the manifold at a fixed height does not 
work, for dimensional reasons. Our pinning method is illustrated in the (2+l)-dimensional 
case in Fig.^and is also in the simulations extended to (3+l)-dimensions. 

We have shown earlier that the average manifold energy develops a logarithmic dependence 
on L z [171118) due to extremal statistics. Since the width of a directed manifold grows as IS it 
is expected that there is a number N z of quasi-independent valleys in the energy landscape [191 
I20U21| . proportional to N z ~ L z /L^ . N z thus presents a scaling parameter, and the Gumbel 
limit is obtained for L fixed in the limit L z — ■* oo, N z — » oo. Hence the logarithmic dependence 
on L z results straightforwardly from extreme statistics. We first study P(E) in the case of 
the "one valley" systems, N z — 1, and then in the case when L z or N z is varied. 

Figs. 13a) and|3b) demonstrate the phenomenology of P(E) for various dimensionalities in 
the "one valley" case. We see that the scaled probability distributions are rather independent 
of D and n though the actual scaling exponents ( and 8 vary greatly. This is most striking if 
one considers the integrated distribution [Fig.^b)]. However, there we concentrate in the part 
of the distribution, which is close to mean, and as the actual shape is usually close to Gaussian 
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Fig. 2 - (a) The normalized energy distributions, (E) = and AE = 1, for manifolds with various 
D and n in a "single valley" case, i.e., pinned when D — 1 or using "cone" geometry when D > 1 
(L z — L). The number of random configurations TV = 5 x 10 and the lateral size of the system 
L = 100 for (D + n) = (1 + 1). Similarly TV = 2.3 x 10 5 and L 2 = 75 2 for (D + n) = (2 + 1); 
TV = 3.2 x 10 5 and L 3 = 25 3 for (D + n) = (3 + 1); TV = 2.5 x 10 6 and L = 70 for (D + n) = (1 + 2); 
TV = 5.1 x 10 5 and L = 50 for (D + n) = (1 + 3). (b) Cumulative distribution / P(E)dE for the 
same data as in (a). The inset shows skewness, ([(E — (£'))/A£] 3 ), versus lateral size, L x — L y — L, 
for (D + n) = (2 + 1). The "cone" geometry is used in the transverse direction. The solid line is the 
least-squares-fit to the data leading to a value of —0.28 ± 0.01. The number of random configurations 
TV = 10 4 for all the other L except L — 75 for which the data is the same as in (a) and (b). 



the differences in the tails cannot even be visible. One of the features of the behavior of P{E) 
is that while the skewness, i.e., 3rd moment, 03 = ([(E — (E)) / AE] 3 ) , is L-dependent the tail 
exponents are much less sensitive to the system size. The inset demonstrates an extrapolation 
of skewness for (2+l)-dimensions, and we obtain 173 = —0.28 ± 0.01. 

Let us now consider what happens as a function of N z = L z /w. If N z <C 1 the interface is 
confined to a "box" which is smaller than the expected roughness, thus criticality is destroyed 
and P(E) becomes a Gaussian. If N z is increased (by L z ) the outcome is depicted in Fig.|3Ja). 
For all the varying dimensionalities we obtain a steady but very slow cross-over of P(E) to- 
wards a fast decay on the side of large energies (E 3> (E)) and a slower one on the other side. 
In the asymptotic limit one expects that the Gumbel distribution is found, which has the 
form (in scaled units) Paumbei ~ exp (E — eupE). The convergence is an example of penulti- 
mate extremal distributions. For N z finite the extremal probability-density-function (PDF) is 
given by a form which has a point- wise limit in the asymptotic distribution |22ll2."Sj . but may 
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Fig. 3 - (a) Stretching exponents, r\- and r/+, of the left and right-hand tails, respectively, of the 
energy distributions, exp{ — \(E — (E)) / AE\ V± }, for (D + n) = (1 + 1) and (2 + 1) versus the system 
height L z in unpinned (and non-cone-geometry) cases. Open circles and filled diamonds are from 
(2 + 1) dimensional data. The first one has N = 2 x 10 5 , while the latter has N = 5 — 8 x 10 4 . For 
both data set L x = L y = L — 40. (D + n) — (1 + 1) is plotted as a dashed line. In that case L = 100. 
In order to show the similarity of the trends with (D + n) — (2 + 1) L z for (1 + 1) are scaled by a 
factor of 0.45. The solid line shows the stretching exponent for the distribution of the minimum of 
several Gaussian distributions min{Xk\X ~ Gauss, k — 1, 2, ... , No} versus the number of Gaussian 
distributions No- In that case the i-axis is scaled so that No = 5.1L Z . The horizontal lines show the 
stretching exponent values rj± of the distributions plotted in Fig.|5Ja)- (b) The skewness values of the 
same distributions as in (a). For (D + n) — (2 + 1) only the data with better statistics N — 2 x 10 s 
is shown. For (D + n) — (1 + 1) and for the Gaussian case we use the same x-axis scalings as in (a). 
The horizontal lines show the skewness values of the distributions plotted in Fig. |5J a ) • ( c ) The scaled 
roughness w/L^ vs. scaled system height Lz/Lr , where £ = 0.41, for (D + n) = (2 + 1) and various 
lateral lengths of the system L x — L y — L — 10, 20, 40, 60, and 80. The number of realizations varies 
from N = 1000 for L = 10 to TV = 600 for L = 80. 
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Table I - The skewness aa and stretching exponent values r\- and n+ of the distributions plotted in 
Fig. ffla). 



D 


n 


0-3 


V- 


V+ 


3 


1 


-0.22 ±0.01 


1.73 ±0.10 


2.25 ± 0.10 


2 


1 


-0.24 ±0.01 


1.79 ±0.10 


2.29 ± 0.10 


1 


1 


-0.26 ±0.01 


1.72 ±0.10 


2.33 ±0.10 


1 


2 


-0.38 ±0.01 


1.58 ±0.10 


2.55 ±0.10 


1 


3 


-0.41 ±0.01 


1.49 ±0.10 


2.60 ±0.10 



e.g. exhibit stretched exponential tails, such that the behavior of the exponents is in agreement 
with the convergence to the asymptotic. Thus ?y_ increases, and rj + decreases. Penultimate 
distributions depend on many factors and can resemble more the other two asymptotic forms 
(than the Gumbel) and the convergence can be very slow. 

Next, one can pose the question: is the one- valley ensemble (with appropriate boundary 
conditions) equivalent to what one gets in the limit N z ~ 1, i.e., are the P's the same? This 
is not a trivial one since the iV 2 -ensemble is "grand-canonical", the number of "valleys" is 
not defined in any sample except in the average sense (recall also the KPZ-growth analogy 
for DP's). Our numerical results imply that this is roughly the case. We also show that the 
convergence of the exponents is reproduced [as depicted in Fig. [3[a)] by a simple Gaussian 
model. We mimick a single- valley distribution by picking the minimum out of a small number 
of independent-and-identically-distributed Gaussian variables. Then N z is varied and the 
smallest of N z such ones is chosen. This gives rise to the last set of data in the figure, and as 
can be seen it reproduces the manifold results qualitatively and almost quantitatively. This 
analogy reinforces the picture of independent valleys and N z as a scaling parameter. 

One may thus consider the tail exponents rj_ and t]+ to be "hyper-universal" as the actual 
values depend mostly on the ensemble and, within the numerical accuracy, almost not at 
all on the dimensionality. The actual values of the tail exponents and the skewnesses are 
reported in Table |U However the actual scaling functions of P(E) are not quite the same, as 
is demonstrated by Fig. |3fb) where the skewnesses of the scaled P(E)'s are shown for "single 
valley" ensembles, and compared to the skewness for various N z or L z . As expected from the 
Gaussian-to-Gumbel cross-over the skewnesses vary continuously. We finish with Fig. Etc), 
which shows reasonably conclusively in (2±l)-dimensions, the average roughness, another 
critical property, is little if at all affected by the ensemble. What the figure demonstrates is 
the insensitivity of the amplitude a(N z ) in w = a(7V 2 )L < ', £ = 0.41 to N z once N z > 1. This 
is in contrast to the average energy ^7] . 

To conclude we have here studied the statistics of the ground state energy of critical 
elastic manifolds in random systems. These are the first numerical calculations of the energy 
distribution, when D > 1. We show that the scaling function of the energy is a continuous 
function of the external control parameter, the system "thickness" N z . Thus the properties 
of the tails of P(E) depend crucially on such details, but not so much on dimensionality or in 
other words the actual scaling exponents of the problem. One should note that the expected 
variation of rj- , r/ + is typically logarithmic in N z using usual convergence rates of penultimate 
distribution extreme statistics. Therefore, such exponents have rather "arbitrary" values. 

There are three open problems that we can outline. First, is it indeed mostly a mathe- 
matical accident that e.g. ?y_ ~ 1.5 . . . 1.7 for all D and n in the single-valley ensembles? The 
existence of such tails tells that all critical manifolds might share some general properties with 
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the KPZ-equation related systems that can be treated on an analytical basis. In particular the 
free energy or large deviation functional could be in such a case non-local. Second, it would be 
interesting (but numerically extremely hard) to look at T > the free energy distribution in 
the general case, except say for the (1 + 1) dimensional directed polymers, which is relatively 
easy and has been studied in the literature Third, one particular example of a "critical 
manifold" is given by random field Ising magnets (RFIM) at the bulk phase transition, since 
here one can map the problem to RB domain walls in (D + 1) dimensions. Our results imply 
that the distribution of the ground state energy should prove to be an interesting quantity, in 
three (3D) and four dimensions. This is so in particular since there is in the much-debated 3D 
case substantial numerical evidence that the aspect ratio of the system will affect the domain 
structure (and thus perhaps the energy distribution). 
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